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The dominant contribution of the next-to-leading order perturbative QCD corrections to the T 
Xy photon energy spectrum for photon energies ft) ^ m[,as can be obtained trivially from the 
knowledge of the corresponding ^7{a) correction to the positronium decay spectrum in the same 



jjj ' region. The latter can be obtained from the NRQED computation of the orthopositronium photon 

spectrum, or by applying the threshold expansion to the 1-loop QED amplitude, as it is explained 
in this talk. 
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The Y — > Xy photon spectrum in the soft-energy region 



Apart from being of physical relevance in itself, the Positronium system provides a testing 
ground for the effective field theory (EFT) concepts and techniques devised in the study of radiative 
decay amplitudes, that could eventually be applied to the description of the more intricate radiative 
decays of heavy quarkonia {QQ —>■ Xy). An example of the latter is provided in this talk, where 
we analyze the region of photon energies ft) ~ ma (soft-energy region) in the radiative decay of a 
heavy fermion-antifermion system through the study of the 3/ decay of the Positronium spin triplet 
ground state (orthopositronium: o-Ps). In particular we show that in this regime the next-to-leading 
order (NLO) corrections to the direct decay photon spectrum of a color singlet QQ pair in angular- 
momentum configuration ^Si are just given by the corresponding correction to the o-Ps photon 
spectrum up to a trivial color factor. The soft-energy region is accessible in the o-Ps^ 3/ decay 
because one of the final state photons can have an arbitrarily small energy, the other two being hard 
photons with energy < m. The process can be then viewed as the radiative version of the o-Ps— > X 
decay, where X = Ij. 

The computation of the iff (a) corrections to the o-Ps spectrum in the soft-energy region can 
be considerably simplified using EFT methods or through the asymptotic expansion of the relevant 
1-loop diagrams. While the latter technique generalizes straightforwardly to the corresponding 
computation for quarkonium radiative decays, the EFT framework naturally explains the origin of 
the dominant term as a consequence of the non-relativistic nature of the system. 

1. NRQED computation of the o-Ps decay spectrum for w ^ m 

The Non-Relativistic QED (NRQED) framework used in Refs. [1,2] provides a systematic 
way to compute bound state effects in the the o-Ps decay spectrum when the photon energy is 
much smaller than the electron mass. Contrary to the usual relativistic approach where the 3- 
photon annihilation is considered to take place at very short distances as compared to the range of 
the electromagnetic binding force between e^e^ , the non-relativistic description takes into account 
that there is a long-distance part in the o-Ps decay process when one the photons in the final state 
is not hard (co <^ m). In the latter case, the decay proceeds in two steps: the low-energy photon is 
first radiated from the bound state, making a transition from the C-odd ground state o-Ps i^Sy) to 
a C-even Positronium state n, which subsequently decays into two photons (see the EFT graph of 
Fig. 1). The emission of the low energy photon from o-Ps is described by the Coulomb Hamiltonian 
of the e^e^ system in interaction with a quantized electromagnetic field [1,2]. 

Since the computation of the photon spectrum in Ref. [1] was intended for photon energies 
comparable to the Positronium binding energy, (O ~ ma^, the interaction Hamiltonian was used 
in the dipole approximation limit. This approximation amounts to evaluating the photon vector 
potential A in the center of mass of the Positronium system, which is fully legitimate for radi- 
ated photons with wavelengths much larger than the characteristic size of the Positronium atom 
{a = 2 /ma). Higher order multipoles arise as a Taylor series in the relative coordinate x and are 
suppressed by powers of ftj/t? under the assumption that the relevant amplitudes one has to compute 
involve integrations to spatial extents of order ~ a. For photons with larger energies, (O 3> ma^, 
this premise will invalidate the use of the multipole expansion. 

However, in the case of the o-Ps system that undergoes a radiative transition before decaying, 
and as it was properly pointed out by Voloshin [3], the scale \/a does not constraint the maximum 
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Figure 1: NRQED graph for the o-Ps^ 3/ decay. The zig-zag Hne represents the low-energy photon. 




Hc + —]G{x,y,K) = 5{x-y), (1.1) 



photon energy for which we can apply the multipole expansion. The reason is that after the soft 
photon is radiated we have to consider all possible e^e^ intermediate states with the right quantum 
numbers. The propagation of these e^e^ states is described by the Green function obtained from 
the Coulomb Hamiltonian, 

m 

at energy K^/m ~ CO —Eq. Since the annihilation of the intermediate e^e^ pair into two photons 
takes place at small distances, the amplitude for the full process is thus given by a convolution 
of the o-Ps ground state and the Green function. In coordinate space, the Green function has a 
characteristic size ruled by the exponential factor exp(— kt). We note that for soft photons of 
energy ma^ <^(0<^m, this exponential factor constrains the relevant integration region to distances 
of order K"^^ ~ \/^/m(0, much smaller than the spatial extent of the initial Positronium atom. 
Therefore, the characteristic distance that enters in the multipole expansion of the soft-photon 
emission in the o-Ps— > 3/ decay amplitude is determined by the falloff of the intermediate Green 
function rather than by the size of the Ps atom, so the series of multipoles is actually an expansion 
in powers of (Or ~ (ajm which can be used as long as ft) ^ m. Note also that for the intermediate 
pair, an iterative computation of the Coulomb Green function, G = + G^VcG^ + . . ., with G^ 
the free Green function, shows that adding a Coulomb exchange generates a term 

f am [m ^ ^ o 

I d^xVcG' r ~aw— <1, for ma^<aXm, 

J K \J (0 

as G^ ~ mexp{—Kr) / r. Therefore the Coulomb interaction can be treated perturbatively in the vir- 
tual e^e^ system after soft photon radiation. The soft-energy region ma^ <^ co <^m thus separates 
the regime where binding effects become essential, from the hard-energy region, where the details 
of the bound state dynamics are irrelevant. 

The NRQED result for the photon energy spectrum in this region, up to NLO reads [2, 3] : 



dPoPs 2ma° 



dx 21n 



14 r- 
5x — — a\Jx + &{ax) 



(1.2) 



with X = (o/m. The LO term of the NRQED specti'um matches the ft) — > limit of the ti"ee-level 
QED calculation first done by Ore and Powell [4]. The a^/x correction also agrees with the leading 
term in the ftj ^ limit of the 1-loop QED result, as shown in Ref. [2] by expanding the analytic 
expression of the exact 1-loop phase-space distribution given by Adkins [5]. 
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Figure 2: Ladder and double- vertex graphs contributing to the o-Ps decay amplitude to order a. 



2. Threshold expansion of the o-Ps^ By 1-loop diagrams 

The NLO term in the NRQED result above arises from distances in the e^e^ system of order 
r ~ 1 / \/mO), that are much smaller than the size of the Ps atom r ~ \/ma. Translated into 
momentum space, this implies that there is a kinematic region where the relative 3-momentum of 
the non-relativistic e^e^ pair obeys the non-relativistic relation ~ p^/m ~ ft). This scaling does 
not correspond to any of the known modes that have been identified in the common applications 
of non-relativistic EFT's, and should be compared to the scaling of the familiar potential modes, 
~ p^/m ~ ma?, characteristic of the heavy particles that form the bound state, and that of 
the soft modes, p° ~ |p| ~ ma » p^/m. The new loop-momentum region has to be taken into 
account for a successful application of the threshold expansion method [6] to QED and QCD loop 
diagrams involving a particle-antiparticle system decaying through the emission of a soft photon. 
In the conventional loop expansion, soft photon radiation from a heavy particle-antiparticle system 
introduces a soft energy component (O ^ ma in the zero component of the heavy particle momenta. 
Non-relativistic poles in massive propagators of the form [p^ — p^/2m — G))^' can be thus found in 
these loops, giving rise to a contribution from the loop-momentum region pQ ^ p^/m ^ ft), which 
we refer to as soft-radiation region. 

To illustrate the latter consider the threshold expansion to the 1 -loop o-Ps^ 3 y amplitude when 
one of the photons has an energy (0\ ~ ma, where the conventional QED perturbative expansion 
is still valid. For our purposes we can take the electron and positron inside the bound state to be 
at rest. Then ftJi represents the smallest scale in the diagrams, and the leading term in the (0\/m 
expansion is given by the ladder and double-vertex graphs shown in Fig. 2. The ladder amplitude 
involves the computation of the 5 -point scalar integral: 

J ^ [ [d^ 

° J p^ [{P+P/2Y - m^\ {{p - P/iy - m2] [{p+P/1 -k,Y- m2] [{p - P /2 + hY - m^\ ' 

with P = (2m, 0) the Positronium momentum. When the loop momentum is hard (p ~ m) the 
corresponding contribution scales as 1^^^ ^ m^^ . However, for ft)i — > the loop-momentum regions 
where p <^ m give larger contributions. Expanding the propagators for p <^ m we have 

^(small) ^ l_ f [dl^ 

^ Im^ J (p2_p2^(2wipo-p^)(-2wipo-p^)(2m7?o-2mft)i -p2) 

which receives contributions from the poles at p^ = (— p2/2m + /£) and po = (— |p| when we 
close the integration contour in the upper complex po-plane. From the residue of the former pole 
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Figure 3: Leading-order Feynman diagrams for radiative T decay, 
we obtain the contribution from the soft-radiation region where /jq ~ p^/m mcoi : 

which gives the leading term in the coi/m expansion of Iq. The pole at po = (— |p| +is) gives an 
integral dominated by p ~ mcoi which scales as [d^p]/{mp())^ ~ m^^co^^, i.e. larger than the hard 
region one, although suppressed by G)\/m with respect the term from the soft-radiation region. 
It is straightforward to verify that only the 1-loop graphs in Fig. 2 receive contributions from the 
soft-radiation region (see Ref. [2] for a more detailed discussion). Once propagator numerators 
and vertex factors are included, the soft-radiation piece in the threshold expansion of the 1-loop 
e^e^ — > 37 amplitude reproduces the NLO term in the photon spectrum shown in Eq. (1.2). 



3. The direct photon spectrum in the radiative T decay for to ~ m^a^ 

The photon energy spectrum in heavy quarkonium radiative decays receives contributions from 
two different mechanisms. The direct photon spectrum entails those contributions where the ob- 
served photon is radiated off a heavy quark, Fig. 3(a), while the fragmentation contributions are 
those in which one of the partons fragments and transfers a fraction of its momentum to the pho- 
ton. Fig. 3(b). The NRQCD factorization approach [7] allows to write down an operator product 
expansion for the direct photon spectrum in, for example, T decays: 

^ = £C„(^)(T|€?„|T), (3.1) 
ax „ 

where the sum above extends over all QQ[n] configurations that can be found inside the quarko- 
nium, and the C„ {x) are short distance Wilson coefficients that can be determined from the anni- 
hilation cross section of the on-shell QQ[n] pair as a perturbative series in as{m},). The ^„ in the 
long-distance matrix element are NRQCD operators which are organized in powers of the relative 
velocity of the heavy quarks. At leading order only the color-singlet operator 0'i(^S\), that cre- 
ates and annihilates a quark-antiquark pair in a color-singlet configuration, contributes. The 
nonperturbative NRQCD matrix element is related to the T wave function at the origin 

(r|^i(35i)|T) = (T|0;a';^_p£p,a>p,|T) =2iV,|./^(0)p. 

The matching coefficient at leading order in a, comes from the tree-level QQ —>■ ggj annihila- 
tion [8] (x = CO /mi, in this case). 



Lo. ^ _ 32 aajQl 



X {l-x)x 2(l-;c)^log(l-x) 2(l-;c)log(l 



■ + - - + 



2 .0 ,.n3 ' ^2 



{2-xY {1-x) 



,(3.2) 
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(a) (6) (p) id) 

Figure 4: A sample of the next-to-leading-order diagrams for the direct contributions to radiative T decay. 

which has identical ;c-dependence to that of the Ore-Powell spectrum for o-Ps [4]. Indeed the 
factorization approach applies trivially to QED bound states, where everything is calculable pertur- 
batively and the bound state wavefunction reduces to the Coulomb one. 

The threshold expansion technique used above for the o-Ps^ 3/ decay amplitude, can be ap- 
plied now to calculate the &{as) correction to the short-distance coefficient C\{x) in the region 
where (O ^ m^av A reduced set of the 1-loop diagrams that have to be considered for the vir- 
tual corrections {bb —>■ ygg) is shown in Fig. 3. The abelian topologies are obviously the same as 
those encountered in the e+e^ 3/ 1-loop graphs. Therefore we know that the leading contribu- 
tion in the ft) / expansion comes from the soft-radiation region in the ladder and double- vertex 
topologies. Figs. 3(a) and 3(b). A quick examination of the pole structure from the propagators in 
non-abelian diagrams reveals that they do not receive contributions from the soft-radiation region, 
and thus give subleading terms in (o/m.},. Likewise, real corrections corresponding to tree-level 
diagrams with one more parton in the final state, bb —>■ Jggg, ygqq are not enhanced for ftj — > 0. 
Therefore we can extract the leading &{(Xs) correction to C\{x) when (O ~ m/^a^ from the QED 
result (1.2) with trivial modifications: 

y3aa^ 
27 j ml 



14 



(3.3) 



From the result (3.3) we obtain a simple expression for the NLO correction to the direct spectrum 
in the soft-energy region for the radiative T decay: 

1 dr^LO 56 1 ^ /4\ 1 



pLO 
^ dir 



dx 



oPs 



pLO 
^ oPs 



dx 



(3.4) 



27(712-9)' 

The factor 4/3 above is just given by the ratio between the colour factors in the ladder and tree-level 
diagrams. The full perturbative ^(tts) corrections to the T photon energy spectrum were calculated 
numerically in Ref. [9]. We have superimposed the NLO result for the spectrum including the 
correction (3.4) to the complete ff{as) result from Ref. [9] in Fig. 3, where = 0.2. As expected, 
our result can only improve the LO prediction for x < 0.3, as it was obtained assuming the scaling 
X ~ a,. The agreement with the exact result for small x is remarkably good considering that our 
formula neglects terms of order ^/x. The approximated NLO direct spectrum becomes negative for 
X < 0.05. However, at those energies (x < aj) the emission of the photon can produce transitions 
to virtual bound states and does not longer belong to the short-distance part of the decay [1]. 
Moreover, we should be aware that the fragmentation contributions to the photon spectrum become 
important in the low-x region [10, 1 1] (namely for x < 0.3 in T decays), and have to be properly 
taken into account for a computation of the full spectrum. The standard NRQCD factorization also 
breaks down at large values of the photon energy [12], where one needs to consider also coUinear 
degrees of freedom [13]. 
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Figure 5: The direct photon energy spectrum in Y radiative decay. The stepped dotted line is the leading- 
order direct photon spectrum, and the solid stepped line includes the full 6'{(Xs) corrections from Ref. [9]. 
The solid blue line is the approximated result given by the leading term in the (a/nii, expansion of the ^(cCj) 
corrections. 
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